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Abstract 



We consider the Maldacena conjecture applied to the near horizon geometry of a Dl- 
brane in the supergravity approximation and present numerical results of a test of the 
conjecture against the boundary field theory calculation using DLCQ. We previously 
calculated the two-point function of the stress-energy tensor on the supergravity side; 
the methods of Gubser, Klebanov, Polyakov, and Witten were used. On the field theory 
side, we derived an explicit expression for the two-point function in terms of data that 
may be extracted from the supersymmetric discrete light cone quantization (SDLCQ) 
calculation at a given harmonic resolution. This yielded a well defined numerical 
algorithm for computing the two-point function. For the supersymmetric Yang-Mills 
theory with 16 supercharges that arises in the Maldacena conjecture, the algorithm is 
perfectly well defined; however, the size of the numerical computation prevented us from 
obtaining a numerical check of the conjecture. We now present numerical results with 
approximately 1000 times as many states as we previously considered. These results 
support the Maldacena conjecture and are within 10 — 15% of the predicted numerical 
results in some regions. Our results are still not sufficient to demonstrate convergence, 
and, therefore, cannot be considered to a numerical proof of the conjecture. We present 
a method for using a "flavor" symmetry to greatly reduce the size of the basis and 
discuss a numerical method that we use which is particularly well suited for this type 
of matrix element calculation. 
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1 Introduction 



The discovery that certain field theories admit concrete realizations as a string theory on a 
particular background has caused a great deal of excitement in recent years However, 
attempts to apply these correspondences to study the details of these theories have only met 
with limited success so far. The problem stems from the fact that our understanding of both 
sides of the correspondence is limited. On the field theory side, most of what we know comes 
from perturbation theory where we assume that the coupling is weak. On the string theory 
side, most of what we know comes from the supergravity approximation where the curvature 
is small. There are no known situations where both approximations are simultaneously valid. 
At the present time, comparisons between the dual gauge/string theories have been restricted 
to either qualitative issues or quantities constrained by symmetry. Any improvement in our 
understanding of field theories beyond perturbation theory or string theories beyond the 
supergravity approximation is, therefore, a welcome development. 

Previously we showed that Supersymmetric Discrete Light Cone Quantization (SDLCQ) 
of field theories ||, £|, |, || can, in principle, be used to make a quantitative comparison with 
the supergravity approximation on the string theory side of the correspondence. We dis- 
cussed this in two space-time dimensions where the SDLCQ approach works particularly 
well; however, it can in principle be extended to more dimensions. 

We will study the field theory/string theory correspondence motivated by considering 
the near- horizon decoupling limit of a Dl-brane in type IIB string theory ||. The gauge 
theory corresponding to this theory is the Yang-Mills theory in two dimensions with 16 
supercharges. Its SDLCQ formulation was recently reported in ||, and recent work has put 
the use of SDLCQ for this class of problems on a stronger footing J7|. This is probably 
the simplest known example of a field theory/string theory correspondence involving a field 
theory in two dimensions with a concrete Lagrangian formulation. 

A convenient quantity that can be computed on both sides of the correspondence is 



the correlation function of gauge invariant operators [10, 11]. We will focus on two-point 



functions of the stress-energy tensor. This turns out to be a very convenient quantity to 
compute for reasons that are discussed in 0. Some aspects of this, as it pertains to a 
consideration of black hole entropy, were recently discussed in |12 |. In the DLCQ literature, 
the spectrum of hadrons is often reported [0] . This would be fine for theories in a confining 
phase. However, we expect the SYM in two dimensions to flow to a non-trivial conformal 
fixed point in the infra-red. The spectrum of states will therefore form a continuum and 
will be cumbersome to handle. On the string theory side, entropy density and the quark 
anti-quark potential are frequently reported. The definition of entropy density requires that 
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we place the field theory in a space-like box which is incommensurate with the light-like box 
of DLCQ. Similarly, a static quark anti-quark configuration does not fit very well inside a 
discretized light-cone geometry. A correlation function of point-like operators does not suffer 
from these problems. 



2 Correlation functions in supergravity 

The correlation function of the stress-energy tensor on the string theory side, with use of 
the supergravity approximation, was presented in p|, and we will only quote the result here. 
The computation is essentially a generalization of [|H], [lT|. The main conclusion on the 



supergravity side was reported recently in [15]. Up to a numerical coefficient of order one, 
which we have suppressed, we found that 

3 

N 2 

(O(x)O(0)> = — ^ . (1) 

This result passes the following important consistency test. The SYM in 2 dimensions with 
16 supercharges have conformal fixed points in both UV and IR with central charges of order 
N% and N c , respectively. Therefore, we expect the two point function of the stress-energy 
tensor to scale like N% / x 4 and N c / x 4 in the deep UV and IR, respectively. According to the 
analysis of || , we expect to deviate from these conformal behaviors and cross over to a regime 
where the supergravity calculation can be trusted. The crossover occurs at x — 
and x = \fN~c/9YM- At these points, the N c scaling of (]!]) and the conformal result match 
in the sense of the correspondence principle O . 



3 Correlation functions in SUSY with 16 Super Charges 

The challenge then is to attempt to reproduce the scaling relation ([I]), fix the numerical 
coefficient, and determine the details of the crossover behavior using SDLCQ. In order to 
actually evaluate the correlation functions, we must resort to numerical analysis. 

The technique of SDLCQ is reviewed in ||, so we will be brief here. The basic idea of 
light-cone quantization is to parameterize space-time using light-cone coordinates x + and 
x~ and to quantize the theory making x + play the role of time. In the discrete light cone 
approach, we require the momentum p_ = p + along the x~ direction to take on discrete 
values in units of p + /K where p + is the conserved total momentum of the system and 
K is an integer commonly referred to as the harmonic resolution 0. One can think of 
this discretization as a consequence of compactifying the x~ coordinate on a circle with 
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a period 2L = 2irK/p + . The advantage of discretizing on the light cone is the fact that 
the dimension of the Hilbert space becomes finite. Therefore, the Hamiltonian is a finite 
dimensional matrix, and its dynamics can be solved explicitly. In SDLCQ one makes the 
DLCQ approximation to the supercharges, and these discrete representations satisfy the 
supersymmetry algebra. Therefore SDLCQ enjoys the improved renormalization properties 
of supersymmetric theories. Of course, to recover the continuum result, we must send K 
to infinity and as luck would have it, we find that SDLCQ usually converges faster than 
the naive DLCQ. Of course, in the process the size of the matrices will grow, making the 
computation harder and harder. 

Let us now return to the problem at hand. We would like to compute a general expression 
of the form F(x~,x + ) = (0(x~, x + )O(0, 0)). In DLCQ, where we fix the total momentum 
in the x~ direction, it is more natural to compute the Fourier transform and express the 
transform in a spectral decomposed form 

F(P_,x + ) = ^(O(P-,x + )O(-P_,0)) =Y,^(0\O(P-)\z)e-^ + {t \ O (-P^0)\0} ■ (2) 

The position-space form of the correlation function is recovered by Fourier transforming with 
respect to P_ = Kit/L. We can continue to Euclidean space by taking r = y/2x + x~ to be 
real. The result for the correlator of the stress-energy tensor was presented in 0, and we 
only quote the results here: 



F(x-,x + ) = J2^(0\T ++ (Km (jp J ^yK 4 (M i V2x^ r ) , (3) 

where Mi is a mass eigenvalue and K±(x) is the modified Bessel function of order 4. In || 
we found that the momentum operator T ++ (x) is given by 

1 



2 / x + \ 2 ill 4 



T ++ (x) = tr 



(<9_X 7 ) 2 + | (iu a d„u a - i{d^u a )u a ) 



I,a=l... 8 (4) 



where X and u are the physical adjoint scalars and fermions respectively, following the 
notation of ||. When discretized, these operators have the mode expansions 



x 1 = 



iJ -'4L 



-I oo 

7= "E \b?An)e- i7Tnx ~ /L + b^(-n)e i7Tnx ~ /L ] . (5) 
/4L n=1 L J 

The matrix element (L/ir)(0\T ++ (K)\i) is independent of L and can be substituted directly 
to give an explicit expression for the two-point function. We see immediately that the 
correlator has the correct small-r behavior, for in that limit, it asymptotes to 
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On the other hand, the contribution to the correlator from strictly massless states is 
given by 



x 



F(x-,x + ) = Y, 



-<0|T++(K)|z>' 2 

7T 



6 



=o KH 2 r A 



(7) 



It is important that this 1/r 4 behavior at large r not be confused with the 1/r 4 behavior 
that we seek at large r. First of all, there is not supposed to be any massless physical bound 
state in this theory, and, secondly, it has the wrong N c dependence. 

Relative to the 1/r 4 behavior at small r, the 1/r 4 behavior at large r that we expect 
is down by a factor of 1/N C . Since we are doing a large- N c calculation, this behavior is 
suppressed. We can only hope to see the transition from the 1/r 4 behavior at small r to the 
region where the correlator behaves like 1/r 5 . 



4 Discrete Symmetries of the Problem. 

In order to calculate the correlation function we use the expression (0). This means that after 
diagonalizing the Hamiltonian P~ one should evaluate the projection of each eigenfunction 
on the specific state T ++ (—K)\0). The fact that we are only interested in states which have 
nonzero value of such projection leads to significant simplifications. 

One can diagonalize any of the eight supercharges Q~. In the continuum limit, the result 
does not depend on the value of a that one chooses, but in DLCQ the situation is a little 
more subtle. As was shown in ||, while the spectrum of {Q^) 2 is the same for all a, the 
wave functions depend on the choice of supercharge. This dependence is an artifact of DLCQ 
and should disappear in the continuum limit. We refer to for the discussion of this issue. 
Here we will just pick one supercharge (for example, Qi). Since the state T ++ (—K)\0) is a 
singlet under R-symmetry acting on the "flavor" index of Q~, the correlator (|2]) does not 
depend on the choice of a even at finite resolution. 

A significant simplification occurs at this stage. Suppose there exists an operator S 
commuting with both P~ and T ++ (— K) and such that S\0) = Sq\0). Then the Hamiltonian 
and S can be diagonalized simultaneously. From now on we assume that the set of states \i) is 
a result of such diagonalization. In this case, only states satisfying the condition S\i) = Sq\i) 
contribute to the sum in (0), and we only need to diagonalize P~ in this sector. So if 
one finds a large enough set of appropriate operators S, then the size of the problem can be 
significantly reduced. By looking at the structure of the state T ++ (— K)\0) one can conclude, 
given arbitrary permutations P and Q of the 8 flavor indices, that any transformation of the 



4 



form 



aljik) ^ f(I)a% I] (k), 
te(fc) -> g(a)l3 a] (k), 



f{I) 

9{<x) 



±1 
±1 



commutes with T ++ (— K), and that the vacuum is an eigenstate of this transformation 
with eigenvalue 1. The requirement for P~ = {Qi ) 2 to be invariant under S imposes some 
restrictions on the permutations. In fact, we will require that Qi be invariant under S, in 
order to guarantee that P~ is invariant. 



The form of the supercharge from || is 

/•CO 

Qa = I [-ftai^aiikjaiih) + ... + (fr/jj - (3j(3j) a p[...}bl{h)a I {k 1 )a J {k 2 ) + ... (9) 
Here the j3i are 8x8 real matrices satisfying {/?/, flj} = 25/ j. We use a special representation 



for these matrices given in JTi ] . 

Let us consider the expression for Q±. The first part of the supercharge (the one which 
does not include (3 matrices) is invariant under (|8|) as long as g(l) = 1 and Q[l] = 1. We 
will consider only such transformations. In order to analyze the symmetries of the f3 terms, 
let us make the following observation. In the representation of (3 matrices we have chosen, 
the expression Bfj = {j3i0j — f3j/3j^j may take only the values ±2 or zero. Moreover, for 
any pair (J, J) there is at most one value of a corresponding to nonzero B. This fact allows 
us to represent B in a compact form. To do so, we introduce a new object \i defined by 



fJ-u 



a . 
-a . 
0. 



Bfj = 2 

Bfj = -2 

Bfj = for all a . 



(10) 



Our choice of j3 matrices then leads to the following expression for fi: 
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(ii) 



We are looking for a subset of transformations (|8|) that satisfy the conditions g(l) = 1 and 
Q[l] = 1 and leave the matrix /x invariant. The latter property means that 



Q[vp[i]p[j)} = g(»u)f(i)f(J)» 



u 



(12) 



Since the subset of transformations that we seek forms a subgroup R of the permutation 
group Sg x S 8 , it is natural to look for the elements of R that square to one. In the case of 
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S$ x Sg, it is known that products of such elements generate the whole group, and, as we 
will show later, the same is true for R. One can construct all Z 2 symmetries satisfying fll~2|) , 
but not all of them are independent. In particular if a and b are two such symmetries then 
aba is also a Z 2 symmetry. By studying different possibilities we have found that there are 
7 independent Z 2 symmetries in the group R, and we have chosen them to be 
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Using Mathematica we explicitly constructed all the symmetries of the type (|8]) satisfying 
(|j~2l). We found that the group of such transformations has 168 elements, and we have shown 
that all of them can be generated from the seven Z 2 symmetries mentioned above. 

In our numerical procedure we use the Z 2 symmetries in the following way. Since all states 
relevant for the correlator are singlets under the symmetry group R, we join our states in 
classes and treat the whole class as a new state. For instance, the simplest nontrivial singlet 
looks like 

|l) = ^Etr(a t (l,/)at(ir-l,J))|0). (13) 
1=1 

This means that if, during the construction of the basis, we encounter the state a'(l, l)a\K— 
1, 1)|0) it will be replaced by the class representative (in this case, by the state |1)). Such 
a procedure significantly decreases the size of the basis, while keeping all the information 
necessary for calculating the correlator. 



5 Numerical Results 

Our numerical results are presented in Figs. 1(a) and 1(b). Figure 1(a) is a log- log plot of 
r 4 times the correlator versus r, so that a 1/r 4 behavior appears as a flat line and a 1/r 5 
behavior gives rise to a line with slope —1. In Fig. 1(b) we plot the log- log derivative, which 
is computed from explicit differentiation inside the sum and amounts to a replacement of 
K,(M ir ) by MiKsiMir). 

Computing this correlator beyond the small-r asymptotics represents a formidable tech- 
nical challenge. In we were able to construct the mass matrix explicitly and compute the 
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spectrum for K = 2, K = 3, and K = 4. Even for these modest values of the harmonic 
resolution, the Hilbert space contained thousands of states. Previously in [fj] we used this 
spectrum and the associated wave function to calculate the correlator beyond the small-r 
region. In the calculation we present here we have made three improvements which have 
allowed us to expand the space by a factor of approximately 1000. The first and most 
straightforward improvement was to rewrite the code in C++, which simply runs faster 
than the Mathematica code and can be exported to faster machines. The second was to use 
the discrete flavor symmetry to reduce the size of the problem at a given resolution. The 
third improvement is a numerical algorithm that replaces the explicit diagonalization with 
an efficient but accurate approximation. 

This numerical algorithm follows from the observation that the contributions to the 
eigenstate sum are weighted by the square of the projection (i\T ++ (—K)\0). The Lanczos 
diagonalization algorithm |15] will naturally generate the states with nonzero projection if 



T ++ (—K)\0) is used as the starting vector. Let \ui) be the normalized vector proportional to 
T ++ (— K) |0), set b\ = 0, and construct a sequence of normalized vectors \u n ) according to the 
Lanczos iteration 6 n+1 |tt n+1 ) = P~\u n ) — a n \u n ) — 6„|m„-i), with a n = (u n \P~\u n ). The \u n ) 
form an orthonormal basis with respect to which P~ is tridiagonal and easily exponentiated. 
Because all of these vectors are generated by applying powers of P to \ui), only those 
eigenvectors with nonzero projections on \ui) can appear. Although generating a complete 
basis by iteration can yield the exact answer ,[] doing many fewer iterations, even 20, can be 
sufficient to capture the important contributions. Such an approach to the computation of 
a matrix element is related to work by Haydock |16|] and others |17j on matrix elements of 
resolvents. 

Before discussing our results we need to address the question of massless states. Our 
SDLCQ calculation of the spectrum of the (8,8) theory saw massless states 0, and we 
argued that they were not normalizable bound states. The argument in that paper was not 
completely correct but the conclusion remains true. We find that in these massless states the 
number of partons in all the contributions is either all even or all odd depending on whether 
the resolution is even or odd. 

We have not, however, removed these unphysical states from the data sets but rather 
used them to obtain an estimate of the region in r where the calculation breaks down. This 
region is where the unphysical massless states dominate the correlator sum. Unfortunately, 
this is also the region where we expect the true large-r behavior to dominate the correlator, if 
only the extra states were absent. The correlator is only sensitive to the two-particle content 



1 Both this statement about the complete basis and the previous statement about nonzero projections will 
hold only in exact arithmetic. Round-off errors will eventually destroy these relationships as the Lanczos 
iteration proceeds. 
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(b) 
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Figure 1: (a)the Log-Log plot of the correlation function (T ++ (x)T ++ (0)) N 2 {2n h +n s ) 
v.s. r in the units where gyM^c/ft — 1 for K = 3,4,5 and 6 (b) the log-log derivative with 
respect to r of the correlation function in (a). 



of the wave function, and we see in Fig. 1(b) the characteristic behavior of the massless 
states at large r only at even resolutions. In Fig. 1(a) for even resolution, the region where 
the correlator starts to behave like 1/r 4 is clearly visible. In Fig. 1(b) we see that for even 
resolution the effect of the massless state on the derivative is felt at smaller values of r where 
the even resolution curves start to turn up. We use these smaller values to estimate the value 
of r where the large- N c approximation breaks down. We see that the value increases as we 
increase the resolution, as expected. Another estimate of where this approximation breaks 
down, that gives consistent values, is the set of points where the even and odd resolution 
derivative curves cross. We do not expect these curves to cross on general grounds, based 
on work in 0, where we considered a number of other theories. 

A proof of the Maldacena conjecture would show up in Fig. 1(b) as a set of derivative 
curves that approached and then touched the line at —1 as we increased the resolution. 
Convergence in the resolution, K, would appear as a flattening of the derivative curves at 
— 1 for the highest values of K. 

We see that the derivative curves are approaching —1 as we increase the resolution and 
appear to be within 10 — 15% before the approximation breaks down. There is however no 
indication of convergence yet; therefore, we cannot claim a numerical proof of the Maldacena 
conjecture. 
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6 Conclusion 



In this article, we used the SDLCQ prescription for computing the correlation function of 
the stress-energy tensor T ++ , which may be readily compared with predictions provided 
by a supergravity analysis following the conjecture of Maldacena JI|. Such a comparison 
requires non-perturbative methods on the field theory side, and the SDLCQ approach is the 
only numerical method suited to this task. At the present time the calculation gives results 
that are within 10 — 15% of the predicted value; however, higher resolution calculations are 
needed to prove convergence. The results we present here increase the number of states by 
a factor of 1000 relative to || . There are currently available methods that we believe could 
give us another factor of 100-1000; however, we have noted in our analysis of our numerical 
results that most of the contributions to the matrix element come from a very small number 
of eigenf unctions. An analytical understanding of this phenomenon could greatly accelerate 
the calculation. 

Finally, we note that, in principle, we could study the proper 1/r behavior at large r by 
computing the 1/N C corrections. In the past we have computed such corrections in some 
theories. However, in the present case such a computation seems to be a very large project 
indeed. 
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